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ABSTRACT 


Theoretical  studies*  of  plasma  acceleration  by  an  rf  field  gradient  indicated 
that  a  plasma  will  be  accelerated  toward  the  low-field  region  if  the  applied  frequen¬ 
cy  is  greater  than  the  critical  frequency  cf  the  plasma,  and  toward  the  high-field 
region  if  the  frequency  is  less  than  the  critical  frequency.  Calculations  for  the 
acceleration  of  plasma  spheres  and  ellipsoids  of  arbitrary  orientation  with  respect 
to  the  rf  field  have  been  completed.  Experiments  with  plasmas  of  various  densities 
have  demonstrated  both  the  acceleration  and  deceleration  of  the  plasmas  as  pre¬ 
dicted  by  theory.  Mercury  plasma  has  been  accelerated  with  a  140-Mc  rf  field 

5 

gradient  to  a  velocity  of  25  x  10  cm/sec,  equivalent  to  a  specific  impulse  of 
2500  seconds. 


ii 


I.  INTRODUCTION 


The  use  of  charged  particles  for  rocket  propulsion  appears  very  attractive 
and  is  under  investigation  in  a  number  of  laboratories  at  the  present  time.  *  The 
primary  advantages  of  such  a  system  over  the  present  chemical  or  mechanical 
system  are  the  high  specific  impulse  generated  and  the  relatively  small  mass  of 
working  fluid  required.  The  investigations  on  charged  particles  cover  work  on 
both  ion  and  plasma  acceleration.  Ions  provide  thrust  when  accelerated  by  elec¬ 
trostatic  fields.  Several  mechanisms  exist  for  obta4ning  thrust  from  plasmas 

g 

(a  neutral  mixture  of  ions  and  electrons).  The  plasma  may  be  heated  electrically  , 
and  the  thrust  provided  by  ejecting  the  energetic  particles  from  the  rocket.  The 

4 

plasma  may  be  accelerated  also  by  means  of  crossed  electric  and  magnetic  fields 

g 

or  by  rapidly  changing  magnetic  fields .  In  addition,  previous  work  at  RCA  Lab- 

g 

oratories  and  at  other  laboratories  have  indicated  the  possibility  of  both  containment 
13  14 

and  acceleration  *  of  a  plasma  by  a  non-uniform  rf  electric  field.  This  latter 
method  is  especially  attractive  because  it  eliminates  the  need  for  heavy  magnetic 
field  structures  in  the  plasma  case  and  eliminates  the  space-charge  problems 
associated  with  the  diode-type  electrostatic  ton  accelerators. 

The  work  reported  herein  has  been  concerned  with  the  demonstration  of  the  ac¬ 
celeration  of  plasma  by  means  of  an  rf  field  gradient  and  deals  primarily  with 
work  of  the  last  six  months.  During  the  first  sb  months  of  this  study  the  major 
experimental  work  included  setting  up  the  proper  equipment  and  performing  pre¬ 
liminary  experiments.  During  that  period,  a  theory  based  on  a  two-particle 
plasma  model  that  predicted  the  force  exerted  on  an  infinite -slab  plasma  by  an 
rf  field  gradient  was  developed.  During  the  last  six  months,  experimental  work 
has  been  concerned  primarily  with  measuring  plasma  densities  and  with  demon- 


-12 


1 


strating  plasma  acceleration  and  decel  j  rati  on.  The  theoretical  work  has  been 
expanded  to  permit  calculations  of  forces  on  plasmas  with  more  realistic  shapes , 
to  explain  some  of  the  experimental  data  obtained,  and  to  suggest  future  experi¬ 
mental  work. 
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II.  THEORY 


A.  INTRODUCTION 

A  theory  of  plasma  acceleration  by  rf  field  gradient  based  on  energy  consid¬ 
erations  has  been  developed.  In  its  simplest  form,  the  theory  yields  the  same 
acceleration  for  a  plasma  slab  as  the  previous  theory  of  T.  V\  Johnston. 1  The 
theory  is  based  on  the  macroscopic  plasma  p  operties,  rather  than  the  motion 
of  individual  particles  considered  in  Johi  iton’ s  theory.  Fewer  and  more 
straightforward  assumptions  are  made,  and  .he  results  have  h-en  shown  to  be 
valid  for  a  wider  range  o~  conditions .  The  theory  can  easily  be  extended  to 
other  geometries;  the  acceleration  of  plasma  ellipsoids  has  been  calculated, 
and  the  results  are  given.  The  theory  can  also  be  extended  to  non-uniform 
plasmas,  and  the  possibilities  and  advantages  of  such  extensions  are  discussed 
in  section  C.  The  magnitude  of  the  effect  of  certain  assumptions  on  tne  theory 
has  also  been  considered.  The  effect  of  a  dc  bias  in  addition  to  the  rf  voltage 

on  the  electrodes  has  been  considered  theoretically  and  the  results  are  presented 
in  Appendix  A. 

B.  ACCELERATION  OF  A  PLASMA  SUB 

The  velocities  attainable  with  rf  field  gradient  acceleration  can  be  calculated 
on  the  basis  of  energy  considerations  as  follows:  The  energy  density  of  an  rf 
electric  field  of  amplitude  E  is 


where  tQ  is  the  permittivity  of  free  space.  When  a  plasma  in  the  form  of  a  thin 
slab  with  boundaries  perpendicular  to  the  electric  field  is  introduced  into  the  field, 
so  that  the  displacement  D  remains  constant,  the  electric  energy  density  inside 
the  plasma  is 


(2) 


where  E ,  as  defined  above,  i^  the  field  present  before  the  plasma  was  introduced. 
The  magnetic  energy  need  not  be  considered  if  there  is  no  applied  magnetic  field 
and  if  the  rf  field  can  be  considered  a  quasi-static  electric  field.  The  dielectric 
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constant  K  of  the  plasma  is 


K  -  —L 


2  2  . 

(O  -  O)  +  I  <ti  V 


-<U  +  *  to  V 


where  <uio  the  frequency  of  the  rf  field,  v  the  electron  collision  frequency,  and 
atp  the  plasma  frequency.  The  plasma  frequency  is  defined  by 


The  potential  energy  U  of  the  plasma  ia  equal  to  the  difference  between  the  electric 
energy  density  inside  the  plasma  and  the  energy  density  present  without  the  plasma: 


U 


P  P 

(6)2-6)2p+((ovp 


(5) 


It  is  assumed  in  what  follows  that  the  plasma  density  is  constant,  that  the  shape 
of  the  plasma  slab  is  not  determined  by  the  field,  and  that  the  motion  of  the  par¬ 
ticles  in  the  plasma  is  small  during  each  period  of  the  rf  field.  If  the  plasma 
moves  to  a  zero  field  region  and  if  the  change  in  electric  energy  density  is  trans¬ 
formed  into  translational  energy,  then 


V 


1  2  12 
—  P  vr  -  -i  p 

2  2 
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where  p  is  the  mass  density  of  the  plasma,  vy  is  the  final  velocity  and  vi  is  the 
initial  velocity.  Thus  the  plasma  attains  a  final  velocity  at  zero  field  of 


As  defined  above,  E  is  the  applied  field  and  therefore  is  independent  of  changes 
in  the  actual  electric  field  due  to  the  plasma.  It  should  be  noted  that  this  final 
velocity  will  be  greater  than  the  initial  velocity  only  if  the  applied  frequency  is 
greater  than  the  plasma  frequency. 


For  comparison  with  theories  of  plasma  confinement,  the  mechanical  force  exerted 
on  the  plasma  can  be  obtained  by  calculating  the  variation  of  the  potential  energy 
V  with  respect  to  a  small  displacement.  The  force  (f)  on  the  plasma  is 


f 


-Vt/ 


(a>*-a>2  )2+(<m>)2 


(8) 


The  mechanism  of  the  energy  transfer  can  be  understood  physically  by  considering 
the  polarization  of  the  plasma  in  the  electric  field  gradient.  For  high  density 
plasmas  (oip>o>)  ,  the  polarization  P  is  in  the  same  direction  as  the  applied  field, 
and  therefore  a  force  P-  Vff  exists  in  the  direction  of  the  field  gradient.  For 
low  density  plasmas  (o>p<  «u)  ;iie  polarization  is  in  the  opposite  direction  and  the 
force  is  reversed. 


C.  EXTENSIONS  OF  THE  THEORY 

One  of  the  principal  advantages  of  the  theory  based  on  energy  considerations 
is  the  possibility  for  extensions  to  more  general  cases.  The  acceleration  of  a 
plasma  of  any  geometry  can  be  calculated  by  substituting  a  dielectric  for  the 
plasma  and  calculating  the  static  electric  field  configuration.  An  extension  of 
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the  theory  to  plasma  ellipsoids  is  given  in  Appendix  £.  Special  cases  can  then  be 
computed,  including  the  acceleration  of  a  slab,  a  sphere,  and  a  flat  disk  by  varying 
the  ratio  of  major  to  minor  axis.  These  calculations  show  that  the  major  effect  of 
geometry  changes  is  a  change  in  resonant  frequency.  The  resonant  frequency  for 
all  other  shapes  is  always  less  than  that  for  a  slab  perpendicular  to  the  electric 
field.  For  a  long  needle  of  plasma,  parallel  to  the  rf  field,  the  resonant  frequency 
decreases  rapidly  as  the  length-to-width  ratio  increases. 

The  assumption  that  the  plasma  density  is  uniform  is  a  difficult  condition  to  realize 
in  practice.  In  the  present  experiment  the  density  varies  by  a  factor  of  100  in  the 
accelerating  region.  In  the  theory  based  on  energy  considerations,  the  problem 
is  equivalent  to  an  electric  potential  problem  with  variable  dielectric  constant. 

While  such  problems  are  difficult  to  handle,  they  can  be  solved  in  principle.  The 
problem  of  concentric  spherical  shells  of  different  densities  has  been  considered. 

One  of  the  principal  results  is  that  resonance  no  longer  occurs  at  a  single  frequency, 
but  occurs  at  several  frequencies. 

D.  ASSUMPTIONS  IN  THE  THEORY 

The  major  theoretical  assumption  is  that  C*e  energy  within  the  plasma  in  the 
rf  field  will  be  converted  into  translational  energy.  For  the  neutral  and  ion  densi¬ 
ties  reported  in  the  following  experiment,  the  mean  free  path  is  too  long  for  scatter¬ 
ing  or  turbulence  to  occur.  Hence  the  field  energy  is  not  lost  as  heat.  Another 
possibility  is  that  the  field  energy  is  converted  to  an  alternating  motion  of  the  elec¬ 
trons  rather  than  an  overall  drift.  The  fraction  of  the  field  energy  converted  to 
drift  depends  on  how  the  rf  field  varies  from  the  maximum  field  region  to  the  min¬ 
imum  field  region.  If  the  field  drops  abruptly ,  the  energy  will  not  be  converted 
into  translational  energy.  In  the  latter  case  the  plasma  oscillations  will  continue 
in  the  zero  fitld  region  and  the  energy  will  eventually  be  converted  into  thermal 
energy.  The  critical  field  gradient  is  not  determined  by  the  present  theory,  but 
In  Appendix  C  it  is  shown  that  the  limit  is  approximate  • 
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To  accelerate  a  plasma  which  starts  in  a  zero  field  region,  the  field  gradient 
should  be  greater  than  the  critical  value  until  the  point  of  maximum  field  is 
reached,  so  that  the  plasma  is  not  decelerated  by  the  rf  field;  then  the  rf  field 
should  decrease  at  such  a  rate  that  the  field  gradient  is  less  than  the  critical 
value,  so  that  the  plasma  can  be  effectively  accelerated.  In  Appendix  C  the 
effects  on  plasma  acceleration  as  the  field  gradient  approaches  the  critical 
value  are  investigated. 

E.  CONCLUSIONS 


Plasma  velocities  that  can  be  attained  by  rf  acceleration  are  given  by 
Equation  7.  The  final  theoretical  velocity  for  a  mercury  plasma  can  be  plotted 
as  a  function  of  the  applied  frequency,  the  result  being  the  resonance  curve 
shown  in  Figure  i.  The  effect  of  the  initial  velocity  is  to  increase  the  final 
velocity,  and  for  simplicity  an  initial  velocity  of  zero  was  chosen.  The  effect 
of  the  electron  collision  frequency  is  a  decrease  in  the  maximum  velocities  near 
resonance,  but  little  change  at  frequencies  for  from  resonance,  as  demonstrated 
fcy  the  two  curves  in  the  Figure.  Calculated  electron  collision  frequencies  for 
the  experimental  plasmas  are  much  lower  than  the  values  chosen  for  illustration, 
but  it  is  expected  that  density  variations  not  considered  in  the  theory  will  reduce 
the  maximum  resonance  and  effectively  increase  the  total  damping.  In  terms 
of  the  present  theory  the  maximum  final  velocity  occurs  at  v 

and  has  the  value 


f(max) 


Vj2  + 


2v2 


t*E 


2  m  it  ( 2<OpV  +  v2) 


1/2 


(10) 


7 


20-r 


RATIO  OF  APPLIED  FREQUENCY  TO  PLASMA  FREQUENCY 

(w/wp) 

Figure  1.  Theoretical  Final  Velocities  of  Plasma  as  a  Function 
of  Applied  Frequency 
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At  this  point  the  acceleration  is  toward  the  field  minimum,  and  this  is  the  optimum 
condition  for  practical  applications. 

For  Figure  1  the  magnitude  of  the  applied  rf  field  was  chosen  as  200  v/cm,  a  value 
easily  attainable.  The  dependence  of  final  velocity  on  the  magnitude  of  the  applied 
nf  field  is  illustrated  in  Figure  2.  If  the  initial  velocity  is  zero  the  final  velocity 
is  proportional  to  the  electric  field.  When  the  initial  velocity  is  appreciable,  the 
curve  is  a  hyperbola,  the  slope  of  the  asymptote  being  a  function  of  <•>  /<uc  as  Bhown 
in  the  Figure.  In  Section  V  the  theoretical  curves  will  be  compared  with  the  results 
Obtained  experimentally. 
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PLASMA  VELOCITY  (10®  cm/sec) 


30 


Figure  2.  Theoretical  Final  Velocities  of  Plasma  as  a  Function  of  rf  Field 
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III.  APPARATUS 


A.  INTRODUCTION 

The  apparatus  described  In  this  report  was  built  to  show  experimentally  that  a 
non-uniform  rf  field  can  exert  a  force  on  a  plasma.  The  apparatus  as  shown  in  Fig¬ 
ure  3  consists  of  a  plasma  source  in  an  evacuated  envelope  fitted  with  pls9ma  detec¬ 
tors  and  an  rf  electrode  system. 

B.  PLASMA  SOURCE 

The  plasma  source  supplies  pulses  of  mercury  plasma.  The  pulse  length  is 
approximately  100  microseconds.  The  igniter  type  source,  which  was  developed 
during  the  first  six  months  of  this  study,  consists  of  a  silicon  carbide  igniter  dipped 
into  a  pool  of  mercury.  The  igniter  is  energized  by  discharging  6  joules  from  an 
energy  storage  capacitor.  The  discharge  is  controlled  by  a  thyratron  switch  that 
closes  at  a  specific  voltage. 

The  density  of  the  plasma  pulse  in  the  rf  interaction  region  is  varied  by  adjusting  the 
size  of  the  source  orifice,  A  source  aperture  of  2. 5  mm  diameter  provides  a  max¬ 
imum  density  of  about  1010  ions  per  cubic  centimeter  at  an  axial  distance  of  1  centi¬ 
meter  from  the  source. 

A  source  that  incorporates  a  plasma  shutter  was  constructed  to  decrease  the  length 
of  the  plasma  pulse.  Such  an  effect  should  provide  means  for  increasing  the  accu¬ 
racy  and  reliability  of  time-of-flight  measurements .  A  diagram  of  the  shutter  is 
given  in  Figure  4.  This  shutter  is  being  tested  at  the  present  time. 
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Figure  3.  Experimental  Apparatus 


L2 


APERTURE 

IGNITER 


GRID#  2  LEAD 


Figure  4.  Electronic  Shutter  for  Plasma  Source 
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C.  THE  RF  SYSTEM 


The  rf  acceleration  system  is  shown  in  Figure  3.  It  consists  of  aplasmasource 
on  the  end  of  the  inner  conductor  of  the  rf  feed  and  looks  into  a  cut-off  waveguide.  The 
electrode  structure  was  designed  to  produce  a  field  maximum  very  near  the  orifice. 

A  measurement  of  electric  field  along  the  axis  both  outside  the  source  and  inside  the 
orifice  was  made  by  the  use  of  a  resistance-board  analog.  The  axial  electric  field 
distribution  as  determined  by  this  method  is  given  in  Figure  5.  The  axial  field  de¬ 
creases  to  1/10  maximum  in  less  than  3  centimeters. 

The  magnitude  of  the  electric  field  E  off  the  axis  is  shown  in  Figure  6.  Because  the 

2 

energy  of  the  plasma  depends  upon  the  square  of  the  electric  field,  cor  tours  of  E  , 
rather  than  E,  near  the  orifice  have  been  plotted. 

Details  of  the  rf  system  are  contained  ' 1  a  previous  report.  Both  a  140-Mc  and  a 
300-Mc  generator  have  been  used  for  the  acceleration  experiments. 

D.  EVACUATED  ENVELOPE  AND  VACUUM  SYSTEM 

The  evacuated  tube  envelope  and  vacuum  system  were  described  in  detail  in  a 
previous  report. 1  Experiments  have  been  conducted  with  the  evacuated  envelope 
continually  pumped.  *‘he  upper  half  of  the  evacuated  envelope  is  covered  with  dry 
ice  to  reduce  the  ambient  pressure  caused  by  the  mercury  from  the  source.  This 
procedure  was  followed  in  all  experiments  to  prevent  an  rf  discharge  in  the  enclo¬ 
sure  during  measurements.  Pressure  in  the  system  is  measured  by  a  Pirani  gauge 
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AXIAL  FIELD  IN  VOLTS  /CM  FOR  I  VOLT  BETWEEN  ELECTRODES 


DISTANCE  FROM  PLASMA  SOURCE  ORIFICE  (cm) 

Figure  5.  Magnitude  of  RF  Electric  Field  On  Axis 
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DISTANCE  FROM  PLASMA  SOURCE  ORIFICE  (cm) 
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and  an  ionization  gauge.  The  vacuum  system  is  capable  of  pumping  down  to  a  pres- 
-6 

sure  of  10  mm  Hg  with  the  system  unbaked. 

E.  DIAGNOSTICS 

The  plasma  source  described  in  this  report  was  designed  to  provide  pulses  of 
plasm  4  less  than  a  millisecond  in  duration.  The  average  and  final  velocities  of  these 
pulses  have  been  determined  by  measuring  their  transit  times  from  the  source  to  a 
probe  and  between  two  probes,  respectively.  Both  single  and  double  wire  probes 
have  been  used.  Each  probe  is  mounted  normal  to  the  axis  of  the  envelope  as  shown 
in  Figure  3. 
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IV.  EXPERIMENTAL  RESULTS 


A.  PLASMA  CHARACTERISTICS 

Several  characteristics  of  the  plasma  pulse  were  investigated  under  conditions 
of  zero  rf  field.  Measurement  of  the  plasma  diffusion  velocity  was  made  previously. 
In  the  last  six  months ,  density  and  electron  temperature  of  the  plasma  were  deter¬ 
mined  in  the  region  near  the  source  orifice.  These  parameters  were  determined 
by  single -probe  measurements.  Various  potentials  with  respect  to  ground  were 
applied  to  the  probe  and  the  peak  of  the  current  was  recorded  as  the  plasma  passed 
by  the  probe.  Plasma  density  and  electron  temperature  were  determined  by  an 
analysis  of  the  resulting  current  voltage  curve.  A  typical  curve  is  shown  in 
Figure  7. 

Electron  temperature  is  determined  from  the  slope,  and  the  charge  density  from 
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the  value  of  the  current  at  which  the  curve  departs  from  an  exponential  dependence 
on  voltage.  Measured  density  contours  near  the  source  with  an  0. 8-mm  orifice  are 
shown  in  Figure  8.  In  Figure  9  are  plotted  curves  of  density  vs.  distance  along  the 
axis  of  the  tube  for  ,  8-mm  and  2. 5-ram  source  apertures. 

The  electron  temperature  of  the  plasma  was  about  20, 000*K  for  a  300-volt  discharge 
in  the  igniter  source.  For  a  discharge  voltage  of  400  volts,  the  electron  temperature 
was  29,000*K,  but  the  plasma  density  remained  substantially  the  same. 

The  probe  system  described  in  Section  in  was  used  to  measure  the  drift  velocity  of 
the  plasma.  Figure  10  shows  the  ion  and  electron  currents  to  the  first  probe  as  a 
function  of  time.  The  electrons  arrive  at  the  probes  before  the  ions,  indicating 
ambipolar  diffusion.  Plasma  velocity  was  determined  by  measuring  the  time  of 
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PROBE  CURRENT  (mo) 


Figure  7.  Single  Probe  Characteristic  I-V  Curve 
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Figure  8.  Plasma  Density  Contour  near  0. 8  mm  Orifice 
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Figure  9.  Plasma  Density  on  Axis  for  Various  Apertures 
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DENSITY  (ions/cm' 


Figure  10.  Oscillogram  of  Ambipolar  Diffusion  Effect 

Ion  and  electron  current  to  1st  probe  as  a  function  of  time.  Lower  curve  is  electron 
current. 

Horizontal  scale:  40 p  sec/div. 

Vertical  scale:  electrons  -  11  ma/div. ,  ions  -  36 p  amp  /div. 

Switching  Frequency  between  ion  and  electron  measurements  is  100  KC 


Figure  11.  Oscillogram  of  Plasi  la  Velocity  (zero  rf) 

Electron  current  to  1st  and  fourth  probe  as  a  function  of  time.  Lower  curve  is 
current  from  1st  probe. 

Vertical  scale:  1st  probe  -  8.  3  ma/div. ,  2nd  probe  -  180 ji  amps/ div. 
Horizontal  scale:  20psec/div. 


flight  of  the  plasma  between  two  probes,  9  cm  apart.  The  oscillogram  in  Figure 

11  shows  the  electron  current  to  the  two  probes  as  a  function  of  time.  The  meas- 

5 

ured  plasma  ambipolar  diffusion  velocity  was  5x10  cm/sec. 

The  ion  temperature  7’+can  be  calculated  from  the  electron  temperature  and  the 
plasma  velocity.  From  the  theory  of  ambipolar  diffusion  as  given  in  Appendix  D, 
the  ion  temperature  is 


in  which  r“is  the  electron  temperature,  v  the  ambipolar  velocity,  M  the  ion  mass, 
and  k  the  Boltzmann  constant.  For  the  measured  values:  T~  =  20,000*K,  v  = 
5xi05  cm/sec,  the  ion  temperature  is  calculated  to  be  440*K. 


B.  APPLICATION  OF  140 -MC  RF  FIELD 

1.  LARGE  SOURCE  ORIFICE  (2.5  MM  DIAMETER) 

As  reported  in  the  previous  semi-aimual  report,  1  probes  were  used  to 
measure  the  plasma  velocity  with  a  140-Mc  rf  field  applied.  The  primary  effect 
observed  was  a  delay  of  the  main  plasma  pulse,  which  arrived  at  the  probe  37 
microseconds  aftei  the  peak  discharge  .Jn  the  absence  of  the  rf ,  and  about  53 
microseconds  after  the  discharge  with  400  volts  rf  applied*  With  the  2. 5-mm 
diameter  orifice  the  electron  current  from  the  plasma  was  observed  at  a  probe 
15.5  cm  from  the  source.  The  plasma  pulse  delay  can  be  seen  in  Figure  i2,  for 
three  different  applied  rf  potentials.  This  last  result  demonstrated  that  the  rf 
field  increases  the  transit  time  to  the  probe,  i.  e. ,  the  rf  field  decelerates  the 

*  A  phototube  aimed  at  the  plasma  source  had  previously  determined  that  the  plasma 
peak  discharge  occurs  about  10  microseconds  after  the  scope  is  triggered  by  the 
breakdown  current. 
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Figure  12.  Oscillogram  of  Plasma  Deceleration 

Electron  current  to  2nd  probe  as  a  function  of  time  with  (a)  zero  RF  potential, 
(b)  200  volts  RF  potential,  and  (c)  400  volts  RF  potential. 

Horizontal  scale:  10 p  sec  div. 

Vertical  scale:  1. 1  ma  div. 
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plasma  emerging  from  the  orifice. 


The  oscillograms  in  Figure  12  indicate  that  there  is  electron  current  to  the  probe 
occurring  within  one  microsecond  of  the  initial  discharge  and  lasting  for  20  micro¬ 
seconds.  This  initial  electron  current  increases  with  increasing  rf  field  and  could 
not  bo  prevented  even  when  the  probe  was  biased  90  volts  negative.  This  current 
is  presumed  to  consist  of  electrons  extracted  from  the  plasma  and  accelerated  down 
the  tube  by  the  rf  field. 

Monitoring  of  the  rf  field  with  a  crystal  detector  while  the  plasma  was  in  the  rf 
interaction  region  showed  about  25%  reduction  in  field  strength  due  to  loading  of 
the  rf  field  by  the  plasma.  However,  the  field  strength  was  sufficiently  high  to 
cause  deceleration. 

2.  SMALL-SOURCE  ORIFICE  (0.8  MM  DIAMETER) 

The  orifice  diameter  was  decreased  to  0. 8  mm  to  reduce  the  plasma 
density  in  the  rf  interaction  region.  The  plasma  was  detected  by  two  probes  located 
at  8  cm  and  22  cm  from  the  source  aperture  on  the  tube  axis.  The  far  probe  is  an 
aquadag  surface  painted  on  the  glass. 

Figure  13  shows  the  electron  current  to  the  probes  as  a  function  of  time.  Oscil¬ 
lations  in  each  pulse  mask  the  location  of  the  plasma  peak,  so  multiple  traces 
were  used  to  average  out  these  oscillations  and  permit  a  time  measurement.  The 
finil  plasma  velocity  was  determined  by  measuring  the  plasma  transit  time  between 
probes.  The  measured  velocity  for  a  maximum  rf  field  amplitude  of  76  v/cm  was 
7x10  cm/sec;  the  velocity  increased  to  18x10  cm /sec  when  the  field  amplitude 
was  raised  to  136  v/cm.  The  initial  sharp  peak,  seen  in  each  oscillogram,  is  an 
electron  pulse  which  occurs  almost  simultaneously  at  each  probe. 
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A  decrease  in  transit  time  with  increasing  rt  fiel '«  was  clearly  observed;  however, 
the  time  at  which  the  plasma  reaches  the  near  probe  is  rather  uncertain  because 
of  oscillations  in  the  detected  signal  and  the  large  initial  electron  pulse.  Measure¬ 
ments  of  the  time  at  which  the  ^lasma  reached  the  aquadag  were  less  subject  to 
these  difficulties.  Since  the  rf  field  decreases  to  i/10  of  its  maximum  value  in 
3  cm  and  the  acceleration  distance  is  small  compared  to  the  total  travel  distance 
of  22  cm,  the  average  velocity  of  the  pla3ma  between  source  and  far  probe  is 
approximately  equal  to  the  final  velocity.  The  time  at  which  the  plasma  emerges 
from  the  orifice  is  taken  to  be  the  time  of  the  initial  electron  pulse.  /erage 
velocities  of  the  plasma  from  the  source  to  the  aquadag  surface  are  plotted  as  a 
function  of  maximum  rf  field  amplitude  in  Figure  14.  The  uncertainties  in  the 
measurements  are  indicated  by  the  vertical  lines  in  the  figure.  There  was  no 
measured  decrease  in  rf  field  as  the  plasma  was  being  accelerated. 

Oscillations  of  the  electron  current  to  the  probes  can  be  seen  in  Figure  13.  The 
frequency  of  the  oscillations  is  about  200  kilocycles  and  increases  with  increasing 
rf  field. 

A  measurement  of  the  plasma  amplitude  at  various  rf  fields  showed  a  definite 
minimum  between  0  and  100  volts /cm  maximum  rf  field  amplitude  as  shown  in 
Figure  15.  The  measurements  were  made  with  a  probe  12.5  cm  from  the  orifice 
using  a  positive  bias  of  45  volts.  The  reproduceability  of  the  curve  was  dependent 
on  the  condition  and  level  of  the  mercury  in  the  igniter  source.  Discussion  and 
explanation  of  this  effect  will  be  found  in  section  V. 
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(VOLT/ CM) 


Figure  14.  Experimental  Plasma  Velocities 
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MAXIMUM  CURRENT  TO  PROBE  I  (  /xamps) 


RF  FIELD  (VOLTS /CM) 


RF  POTENTIAL  (VOLTS) 

Figure  15.  Amplitude  Varlatlo  x  of  Plasma  Pulse  with  Low  rf  Fields 
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V.  DISCUSSION 


A.  COMPARISON  OF  EXPERIMENTAL  RESULTS  AND  THEORY 

1.  DECELERATION 

When  the  source  orifice  was  adjusted  to  a  diameter  of  2. 5  mm  and  a  plas¬ 
ma  pulse  was  produced,  deceleration  was  observed  when  the  140-Mc  rf  field  was 
applied.  Density  measurements,  as  described  in  Section  IV,  show  that  the  plasma 
emerges  with  an  approximately  spherical  shape.  According  to  the  theory,  deceler¬ 
ation  will  occur  in  a  140-Mc  field  for  spherically  shaped  plasma  with  a  density  great¬ 
er  that  the  critical  density 

2 

nc  =  7.5  x  left  *ons /cm  ^  ( from  Appendix  B  nc  =  ^  m(Q  ) 

e ^ 

The  density  of  the  plasma  emitted  from  the  2. 5-mm-diameter  orifice  is  greater 

than  the  critical  density  in  the  most  effective  part  of  the  interaction  region.  The 

plasma  density  decreases  to  the  critical  density  2.4  centimeters  from  the  source, 

the  point  where  the  electric  field  has  decreased  to  15  percent  of  maximum.  The 

5  5 

average  velocity  to  the  prob^  changed  from  4.  2  x  10  cm/sec.  at  zero  rf  to  2. 9  x  10 
cm/ sec.  at  340  v/cm  rf  field.  This  observed  deceleration  is  in  qualitative  agree¬ 
ment  with  the  theory.  In  addition,  the  substitution  of  these  values  in  the  theoretical 
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expression  for  final  velocity  <  equation  6)  yields  a  value  for  density  of  5.  7  x  10 
ions  /cm3.  Such  a  calculated  value  is  consistent  with  the  measured  value  at  a  dis¬ 
tance  of  about  8  mm  from  the  source. 

2.  ACCELERATION 

The  acceleration  predicted  by  theory  has  also  been  observed  experimentally. 
According  to  the  theory,  acceleration  will  occur  in  a  140-Mc  field  for  a  spherically 
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8  3 

shaped  plasma  with  the  density  less  than  the  critical  density  of  7.5  x  10  ions/cm  . 
Plasma  densities  less  than  this  critical  value  were  obtained  with  the  source  orifice 
of  0.8  mm  diameter.  Density  measurements  at  the  peak  of  a  plasma  pulse  show 
that  this  critical  density  occurs  at  a  distance  of  about  0.  6  cm  from  the  source.  At 
all  greater  distances,  the  density  is  less  than  the  critical  density.  At  this  distance 
of  0.6  cm,  the  field  has  only  dropped  to  60%  of  the  peak  so  that  a  considerable  por¬ 
tion  of  the  plasma  is  in  an  accelerating  field  and  thus  acceleration  is  expected.  The 
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velocity  was  observed  to  increase  from  5  x  10  cm/sec  at  zero  rf  field  to 
25  x  105  cm  sec  at  an  rf  field  of  170  v  cm. 

3.  PLASMA  VELOCITY  VERSUS  RF  FIELD 

The  velocity  that  was  obtained  experimentally  using  an  orifice  of  0.8  mm 
diameter  is  plotted  against  rf  field  in  Figure  16.  Also  plotted  are  theoretical  curves 
for  values  of  the  operating  frequency  of  1.05  and  1.025  times  the  critical  frequency, 
with  an  assumed  initial  velocity  equal  to  the  measured  value. 

The  observed  velocity  increases  with  field,  in  qualitative  agreement  with  the  theory. 
To  obtain  the  high  velocities  observed  with  the  magnitude  of  field  used,  theory  re¬ 
quires  that  the  critical  frequency  for  the  plasma  be  very  close  to  the  applied  fre¬ 
quency.  The  fact  that  best  agreement  exists  for  a  value  of  wc/<u  so  near  to  1  sug¬ 
gests  that  either  an  internal  mechanism  exists  for  assisting  the  resonance  or  that 
only  that  portion  of  the  plasma  with  the  critical  density  is  accelerated.  In  either 
case,  such  an  effect  may  be  exploited  in  future  experiments  to  shift  from  a  pulsed 
to  a  continuous  type  of  source. 

The  discrepancy  between  the  observed  and  theoretical  slopes  may  be  attributed  to 
theoretical  assumptions,  experimental  errors,  or  a  change  in  plasma  character¬ 
istic  as  the  field  changes.  Theoretical  assumptions  that  the  plasma  originates  at 
the  rf  field  maximum  and  that  the  plasma  density  is  uniform  were  not  realized  in 
practice.  A  change  in  plasma  characteristic  due  to  focusing  by  the  radial  field 
gradient,  shown  in  Figure  6,  may  increase  the  plasma  density,  and  thus  the 
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Figure  16  Comparison  of  Experimental  and  Theoretical  Plasma  Velocities 
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increased  slope  would  be  expected  because  of  the  decrease  in  the 
ratio  of  <o  /uc  • 

B.  OSCILLATIONS  IN  PROBE  SIGNALS 

For  an  accelerated  plasma,  oscillations  with  a  period  of  about  5  microsecond 
were  observed  in  the  probe  signals.  The  oscillation  frequency  increased  with  in¬ 
creasing  applied  field.  With  the  larger  orifice,  when  the  plasma  was  decelerated, 
this  type  of  oscillation  was  not  observed. 

There  are  several  possible  explanations  for  these  oscillations: 

( 1 )  Plasma  having  a  density  well  below  the  critical  density  will  not  be  appreciably 
accelerated  by  the  rf  field  gradient.  Thus  plasma  emerging  from  the  source  m\*9t 
approach  a  density  close  to  the  critical  density  before  acceleration  takes  place. 

This  plasma  is  accelerated  leaving  behind  a  plasma  void.  More  plasma  moves 
from  the  source  into  the  void  and  the  process  is  repeated.  The  frequency  of  this 
cyclic  process  is  dependent  on  the  plasma  diffusion  velocity. 

(2)  If  theie  is  an  axial  field  gradient  there  must  be  a  radial  field  gradient.  When 
the  plasma  density  is  almost  equal  to  the  critical  density  the  plasma  will  be  strong¬ 
ly  focused  by  the  radial  gradient.  The  focusing  will  produce  a  pinching  effect  so 
that  the  plasma  would  emerge  from  the  rf  interaction  region  as  individual  plasma 
puffs . 

( 3 )  Acceleration  of  the  plasma  results  in  a  leading  negative  space  charge  region 
(electrons)  and  a  trailing  positive  space  charge  region  (ions).  Plasma  moving  in 
behind  the  accelerated  plasma  sees  the  positive  space  charge  region.  The  ions 
are  repelled  and  space  charge  density  fluctuations  may  be  established. 

Oscillations  of  the  above  type  have  also  been  observed  in  experiments  with 
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thermionic  converters  containing  cesium  plasma.  In  the  latter  case,  it 
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has  been  shown  that  these  oscillations  arise  from  an  instability  in  the  plasma 
potential. 

C.  VARIATION  OF  PROBE  SIGNAL  AMPLITUDE 

At  rf  field  levels  below'  120  v/ cm,  the  probe  signal  varies  in  an  interesting 
and  explainable  manner.  As  shown  in  Figure  15,  the  amplitude  of  the  probe  elec¬ 
tron  current  as  a  function  of  rf  field  exhibits  a  minimum  at  45  v/cm. 

To  explain  these  results,  it  is  necessary  to  observe  that  the  plasma  does  not  start 
at  a  point  of  maximum  field,  as  assumed  by  the  theory,  but  starts  inside  the  igniter 
source  at  a  point  of  zero  field.  For  very  low'  rf  voltages  the  plasma  is  decelerated 
until  it  reaches  the  maximum  field  point,  and  then  accelerated.  As  the  rf  voltage 
is  increased,  some  of  the  plasma  will  not  have  enough  initial  kinetic  energy  toover- 
come  this  deceleration,  and  will  be  turned  back  before  it  reaches  the  maximum 
field  point.  For  fields  higher  than  about  45  v/cm,  electrons  can  be  accelerated  all 
the  way  to  the  maximum  field  point  in  less  than  one  rf  cycle  ani  are  not  contained 
by  the  field  barrier.  Hence  the  number  of  electrons  observed  increases  for  higher 
rf  fields.  The  oscillation  amplitude  for  electrons  in  a  140-Mc  field  of  45  v/cm  is 
about  1.0  mm,  which  is  also  the  distance  between  the  maximum  field  point  and  a 
point  in  the  orifice  where  the  field  is  negligible.  Therefore  the  magnitude  of  the 
rf  field  for  which  the  minimum  amplitude  is  observed  is  in  good  agreement  with 
the  theory. 

Appendix  E  contains  a  more  detailed  theoretical  treatment  of  this  effect. 

D  THRUST  CALCULATIONS 

Calculations  show  that  the  thrust  (MNcv2)  from  the  system  at  the  peak  of 
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a  pulse  is  about  0. 13  dynes/ cm  .  Since  the  duty  cycle  is  extremely  small,  the 
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average  thrust  is  too  small  for  measurement  and  a  detector  with  both  sufficiently 
high-frequency  response  and  sensitivity  to  measure  the  transient  thrust  is  not 
available 


For  measurement  purposes,  the  average  thrust  must  be  increased  significantly, 

A  larger  thrust  can  be  obtained  by  increasing  the  plasma  density  and  increasing 
the  driving  frequency  correspondingly.  The  average  thrust  can  also  be  increased 
by  using  a  faster  pulse  rate  or  a  continuous  source  instead  of  the  present  low  duty 
cycle.  The  increased  duty  cycle  in  conjunction  with  a  higher  density  plasma  should 
produce  sufficient  thrust  to  be  detected  by  mechanical  means  and  the  results  of  such 
experiments  would  provide  additional  data  for  comparison  with  the  theory. 

E.  ENERGY  OF  ACCELERATED  IONS 

The  high  translational  energy  achieved  by  the  ions  is  noteworthy.  This  high 

ion  energy  has  been  observed  both  in  the  ambipolar  diffusion  mechanism  and  in 

acceleration  by  an  rf  field  gradient.  In  the  ambipolar  diffusion  experiments,  where 

the  ion  acceleration  is  caused  by  the  pull  of  the  high  speed  thermal  electrons,  the 

5 

ion  velocity  was  5  x  10  cm  /sec.  This  velocity  is  equivalent  to  an  energy  of  25 
volts.  In  the  case  of  the  rf  field  gradient  acceleration,  a  maximum  rf  potential 

5 

of  only  200  volts  was  applied  and  the  velocity  of  the  plasma  was  increased  to  25 x  10 
cm/sec. ,  equivalent  to  an  ion  energy  of  625  volts.  During  the  acceleration  period, 
electrons  must  transfer  the  energy  they  abstract  from  the  field  to  ions  quite  rapidly. 
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VI.  CONCLUSIONS  AND  PLANS 


A.  MAJOR  accomplishment; 

During  the  first  six  months  of  this  study,  theoretical  analyses  predicted 
that  rf  forces  may  be  directed  toward  the  field  maximum  as  well  as  away  from  it 
and  experimental  results  demonstrated  a  change  in  magnitude  of  the  plasma  velocity 
due  to  the  rf  field.  During  the  second  six  months,  the  theory  was  extended  from  a 
two  particle  model  to  one  based  on  macroscopic  plasma  properties.  The  new  theory 
yields  the  same  final  result  as  the  earlier  theory  and  in  addition  permits  calculations 
of  acceleration  for  more  realistic  plasma  shapes  than  the  earlier  assumed  infinite 
slab.  Experimental  results  have  demonstrated  both  the  acceleration  and  deceleration 
of  a  plasma  in  accordance  with  the  theoretical  predictions. 

B.  SOME  DETAILS  OF  THE  ACCOMPLISHMENTS 

Calculations  for  the  acceleration  of  plasma  spheres  and  ellipsoids  with 
any  orientation  with  respect  to  the  rf  field  have  been  completed.  These  show  that 
accelerations  of  plasmas  with  a  density  higher  than  that  in  a  slab  are  possible  for 
a  specific  rf  frequency.  In  particular,  for  a  10-to-l  length-to-width  ratio  in  an 
ellipsoid,  the  density  at  which  the  direction  of  force  changes  is  fifty  times  that  for 
a  plasma  slab.  Hence,  for  a  given  rf  frequency,  use  of  such  a  needle-shaped  plasma 
could  allow  acceleration  of  plasma  of  density  fifty  times  larger  than  is  possible  with 
a  plasma  slab. 

A  reliable  plasma  source  was  developed.  This  source  permits  the 
density  of  the  plasma  entering  the  interaction  region  to  be  varied  by  adjusting  the 
source  aperture.  Plasma  density  contours  were  determined  in  the  interaction  region 
near  the  orifice  for  various  orifice  sizes.  Probes  have  been  used  to  detect  the  moving 
plasmas  and  to  determine  plasma  transit  times.  Plasma  has  been  accelerated  to  a 
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velocity  of  25  x  10  cm/sec  which  i?  equivalent  io  a  speciiic  impulse  of  2500  seconds. 
Low-density  plasnas  were  accelerated  and  high-density  plasmas  were  decelerated 
as  predicted  by  theory.  The  variations  in  plasma  velocity  and  probe  signal  amplitude 
with  changes  in  rf  field  magnitude  were  compared  to  and  explained  by  the  theory. 

C.  FUTURE  WORK 

To  provide  a  quantitative  test  of  the  theory,  the  dependence  of  the  accelera¬ 
tion  on  the  applied  frequency  will  be  determined.  Higher  frequencies  will  permit  the 
use  of  higher  density  plasmas;  this  should  result  in  thrusts  of  sufficient  magnitude  to 
be  detected  experimentally.  Measurements  at  300  Me  are  in  progress.  The  construction 
of  a  2.  5  kMc  system  is  planned.  The  increased  plasma  densities  ior  use  with  this 
high-frequency  system  should  be  available  from  a  source  of  design  similar  to  that 
used  at  lower  frequencies. 

The  igniter  source  will  be  provided  with  an  electronic  shutter  to  define  the  plasma 
pulse  shape  more  sharply.  A  sharply  defined  plasma  pulse  should  assist  in  improv¬ 
ing  the  accuracy  of  velocity  measurements. 

Requirements  of  a  continuous  operation  system  rather  than  a  pulsed  one  will  be  in¬ 
vestigated.  With  experimental  information  from  the  work  at  various  frequencies, 
specific  impulse,  thrust,  density  and  optimum  frequency  will  be  determined. 

Theoretical  studies  of  the  properties  of  plasmas  in  rf  field  gradients  will  be  continued. 
Interpretations  of  the  experimental  results  at  higher  frequencies  and  densities  will 
be  made.  The  theory  for  the  acceleration  of  non-uniform  plasma  is  being  developed 
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APPENDIX  A 


EFFECT  OF  DC  BIAS 

This  section  discusses  the  effect  of  dc  bias  on  rf  acceleration.  It  is  assumed 
that  a  dc  potential  is  applied  on  the  same  electrodes  as  is  the  rf;  hence,  the  rf  and 
dc  field  distributions  are  the  same.  The  analysis  is  similar  to  the  calculations  for 
rf  alone,  *  and  the  same  basic  assumptions  are  made.  A  fully  ionized  plasma  of 
constant  density  in  the  form  of  an  infinite  slab  perpendicular  to  the  rf  field  is  assumed. 
It  is  further  assumed  that  the  plasma  slab  continues  as  a  rigid  plasma  slab.  The 
latter  assumption  may  be  even  less  justified  for  dc  than  it  is  for  rf;  nevertheless, 
the  theory  yields  some  interesting  results. 

The  major  effect  of  a  dc  bias  is  to  change  the  effective  plasma  resonant  frequency. 

The  applied  frequency  at  which  the  acceleration  changes  sign  is  greater  if  the  gradi¬ 
ent  of  the  dc  field  is  positive,  and  less  if  the  gradient  is  negative.  To  decrease  the 
effective  plasma  resonance  the  center  electrode  should  therefore  have  a  positive 
charge. 

Before  proceeding  to  the  theory,  it  may  be  useful  to  visualize  the  process  that  changes 
the  frequency.  Consider  a  plasma  slab  with  no  applied  fields.  The  electron  potential 
energy  can  be  calculated  for  various  positions  of  the  electron  center  of  mass.  If  the 
electrons  are  displaced  from  their  equilibrium  posdions  an  oscillation  will  take 
place,  the  frequency  being  determined  by  the  curvature  of  the  potential  energy  curve 
and  the  mass  of  the  electrons.  If  a  uniform  dc  field  is  apj.  jed,  the  equilibrium 
position  of  the  electrons  is  changed,  corresponding  to  the  polarization  of  the  plasma, 
but  the  resonant  frequency  remains  constant,  as  shown  in  Figure  17.  If  a  dc  field 
gradient  's  applied,  a  curved  potential  is  added,  which  changes  the  curvature  of 


38 


I 

^  UNIFORM  dc  FIELD 

I  V 

I 

ELECTRIC 
POTENTIAL 

z 


I 

j 

i 


ELECTRON-  I  ON  SEPARATIONS  ELECTRON-1  ON  SEPARATIONS 

Figure  17.  Change  in  Plasma  Frequency  bv  D.  C.  Field  Gradient 
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the  potential  energy  curve  and  hence  changes  the  resonant  frequency,  as  shown 
in  figure  17. 


The  differential  e.  tions  of  motion  for  the  plasma  center  of  mass  Zm  and  the 
separation  5  between  the  ion  and  electron  center  of  masses  are*: 


d?Zm 

dt^ 


dE_ 

dz 


(A-i) 


d?8  d8  2S 


e_  dE_  g 
m  dZ 


(A-2) 


dE 

E  is  the  applied  field  and  is  the  applied  field  gradient  at  Zm  ;  v  is  the  elec¬ 

tron  collision  frequency,  the  plasma  frequency,  e  and  m  the  electron  charge 
and  mass,  and  M  is  the  ion  mass. 


The  applied  field  consists  of  an  rf  voltage  and  a  de  bias,  and  can  be  written  in 
the  form 

E  =  Erf  cosail  +  Ejc  (A— 3) 

The  result  of  substituting  equation  (3)  into  equation  (2)  is 
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It  is  convenient  to  define  the  following  quantities: 


CO 


2 

S 


(A-5) 


S' 


(A-6) 


6) 

P 


2 

f 


JL  £jLl 

m  dZ 


(A— 7) 


(A-8) 


The  result  of  substituting  the  above  quantities  into  equation  (A-4)  i3  to  eliminate 
the  constant  term  on  the  right  L.  d  aide: 


— — +  V  +  CO  ?  5  '  +  £J  ‘  5  '  COScjt 

dt 2  dt  P  * 


—  E '  ,  coscot 

m  rf 


(A— 9) 


The  equation  is  now  in  the  same  for»  i  as  the  equation  m  the  theory  without  the  dc 

2 

bias,  and  can  be  solved  by  similar  methods.  Neglecting  the  fourth  terms  (if 

2 

«  oj  ,  )  ,  the  steady  state  solution  of  equation  (A-9)  is 
P 


S' 


_ cos  (ut  -(f)) _ 

If  «  ,2  -  CO~  )~  +  ((Ol>)2} 

P 


% 

(A— 1 0) 


*Note  that  equation  (A-6)  merely  separates  out  the  average  plasma  separation 
which  is  the  polarization  produced  by  the  dc  field.  Similarly  equation  (A-7) 
expresses  the  change  in  plasma  resonance  frequency  due  to  the  dc  field.  The 
latter  equation  is  somewhat  more  symmetrical  if  the  electron  mass  is  not 

compared  to  the  ion  mass,  in  which  case  the  resonance  frequency  is 

given  by: 

(a,  ")  2  =  M’*'  j/  +  m  +  e  (M-m)  r^dc 

P  ( 0  Mm  Mm  SZ 

so  that  the  dc  bias  effect  vanishes  if  the  masses  of  the  positively  charged  particles 
equals  that  of  the  negatively  charged  partirles. 
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where 


tan  <f>  = 


CjJ  l 


(A-U) 


Combining  i  ations  (A-l),  (A-3),  (A-6)  and  (A-8),  the  instantaneous  acceleration  of 
the  i  center  of  mass  is 


it2 


COS cj  t  + 


e  _  L  E.r 

—  75T  rf  - 


COS  — 
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oj- 


2  2\2  /  \2J/2 

+  (coi')  ] 


(A-12) 


The  average  acceleration  is  therefore 


3Z  tmM 


dE 


rf 


a>p'2  -  co? 


dZ 


[(cap  /  2-cj2)2+  (cjv)2  ] 


(A-13) 


The  first  term  is  simply  the  acceleration  towards  the  maximum  of  an  electric  field 
that  any  dielectric  experiences. 


If 


and 


(A-14) 
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<i.'  ,  V 

p 


the  square  of  the  frequency  which  maximizes  (A-15) 
the  acceleration 


up  *  >  >  * 


(A-16) 
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then 


- VErf 

8tnM  \>  (Op* 


(A-17) 


Thus  it  has  been  shown  that  the  average  acceleration  is  in  the  sar\e  form  as  the 
theoretical  result  without  the  dc  bias,  except  for  the  change  in  -  .sma  resonant 
frequency.  If  the  bias  is  increased  too  far,  the  terms  which  were  neglected  would 
become  important,  and  the  acceleration  would  be  towards  the  field  maximum. 
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APPENDIX  B 


ACCELERATION  OF  PLASMA  ELLIPSOIDS 

The  theory  presented  in  the  text  on  the  acceleration  of  a  plasma  olab  in  an 
rf  field  gradient  can  be  extended  to  include  the  acceleration  of  a  plasma  ellipsoid. 
The  theory  will  then  include  the  slab  as  a  special  case,  as  well  as  most  other  cases 
of  interest.  The  energy  of  a  dielectric  ellipsoid  is  calculated,  with  the  dielectric 
constant  that  of  the  plasma.  The  force  on  the  plasma  can  then  be  obtained  by 
differentiating  the  potential  energy  with  respect  to  position,  or  the  final  velocity 
achievable  can  be  calculated  by  equating  the  potential  energy  to  kinetic  energy.  As 
before,  it  is  assumed  that  the  plasma  is  fully  ionized,  of  constant  density,  and  that 
it  is  in  a  quasistatic  electric  field  which  is  approximately  uniform  over  the  distance 
the  plasma  moves  during  one  rf  cycle.  It  is  also  assumed  that  the  shape  of  the 
ellipsoid  is  unaffected  by  the  accelerating  forces. 


The  dielectric  ellipsoid  is  assumed  to  have  semi-principal  axes  of  length  a,  b,  c, 

oriented  along  the  *,  y,  z  axes  respectively,  and  a  dielectric  constant  K  .  The 

1 8 

potential  energy  of  the  ellipsoid  is  then 


V--~  (K-l)t 

h 


OX 


oz 


;+£|£  (K-I)Aj  (K-1)A2  l+2k.(K~l)A. 


(B-l) 


where  E0  is  the  amplitude  of  the  applied  rf  field  and  where 
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l  " 
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ds 


(s  +  a  )R. 
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a2.  f 


ds 


( *  +  b  )  /?, 


(B-2) 


A3-  f - - -  where  Rs  =  f(s+a2)(s  +  b2)(s+c2) 

J0  (s  ♦  c2)  R S 
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The  dielectric  constant  of  plasma  of  constant  density  is  given  by 


15 


K-l  - 


') 


2 

— Cl)  +  /  Cl)  V 


(B-3) 


This  dielectric  constant  is  substituted  into  Equation  B-l,  and  the  real  part  of  V  is 
calculated.  If  this  energy  is  converted  to  translational  energy  the  final  velocity  for 
a  plasma  in  a  field -free  region,  which  had  an  initial  velocity  vi  in  an  electric  field 
E0  is 
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(B-5) 


A  plasma  ellipsoid  with  unequal  axes  has  three  main  resonant  frequencies.  Other 
frequencies  may  be  excited  if  the  applied  electric  field  is  not  uniform.  It  is  useful 
to  noie  a  simple  relation  between  these  frequencies,  namely  that  the  sum  of  the 
squares  is  equal  to  square  of  the  resonant  frequency  for  a  slab,  i.e. , 


2  2  2 
a)  =  a) ,  +  co 
P  1  2 


(B-6) 


This  can  be  proved  by  an  integration  by  parts  of  .4^  .  For  a  plasma  sphere  the 
three  frequencies  are  equal,  therefore  the  resonant  frequency  of  a  sphere  is 


c 


(B-7) 


The  actual  resonant  frequency  in  general  cannot  be  written  in  analytical  form,  since 
the  integration  is  not  possible.  However,  special  cases  can  be  evaluated,  and  a 
number  of  them  are  presented  in  Table  B-i.  Included  is  the  case  of  a  plasma  slab, 
which  has  a  higher  resonant  frequency  than  any  of  the  other  cases.  Of  special  interest 
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TABLE  B-1 


RESONANT  FREQUENCIES  FOR  PLASMA  ELLIPSOIDS 


where 


f  =  b/a  e  =  J  1  -  (b2/a2)  f  =  J~j-(a2/b2) 
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is  the  geometry  with  axial  symmetry,  where  the  electric  field  is  along  the  symmetry 
axis.  The  resonant  frequency  is  plotted  in  Figure  19  as  a  Junction  of  the  ratio  of  the 
third  axis  to  the  other  two. 


For  comparison  to  theories  of  plasma  confinement,  the  mechanical  force  exerted 
on  the  plasma  can  be  obtained  by  calculating  the  variation  of  the  potential  energy 
V  with  respect  to  a  small  displacement.  The  force  on  the  plasma  f_  is  equal  to 


f  =  -  S/U  = 


ne  * 

pn 


(j)  -  co 
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(cj2~  co  ^)2+(cou)2 


2  2 
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oz 


(B-8) 


A  conclusion  from  this  analysis  is  that  acceleration  of  other  plasma  geometries  is 
very  similar  to  that  of  a  plasma  slab;  the  principal  difference  is  in  the  relationship 
between  the  resonant  frequency  and  the  plasma  density.  For  acceleration  to  a  field- 
free  region  it  is  necessary  that  the  applied  frequency  be  above  the  resonant  frequency; 
thus  to  accelerate  a  high-density  plasma  it  is  desirable  to  have  it  assume  a  shape 
such  as  to  minimize  its  resonant  frequency.  The  most  important  conclusion  of  the 
theory  presented  is  that  for  a  long  needle  of  plasma,  parallel  to  the  rf  field,  the 
resonant  frequency  decreases  rapidly  as  the  length  to  width  ratio  increases.  In 
particular,  as  seen  in  Figure  18,  for  a  iO-to-1  length-to-width  ratio,  the  critical 
density  is  fifty  times  that  for  a  plasma  slab.  Use  of  such  a  needle-shaped  plasma 
could  allow  acceleration  of  plasma  of  density  fifty  times  larger  than  is  possible 
with  a  slab  of  plasma,  provided  that  the  assumption  of  a  constant  shape  can  be 
satisfied. 
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BLANK  PAGE 


APPENDIX  C 


EFFECT  OF  URGE  FIELD  GRADIENTS 

Theoretical  calculations  of  the  plasma  acceleration  by  rf  have  been  made  for 
small  field  gradients  In  the  theory  based  on  energy  considerations  the  motion  of 
the  center  of  mass  of  the  plasma  was  assumed  to  be  small  in  each  period  of  the  rf 
field.  In  the  original  mathematical  formulation,  *  a 2  is  assumed  small  compared 
to  «p2  ,  where  o>  2  *  e/m  .  Such  an  approximation  is  useful  in  understanding 
plasma  acceleration. 

The  present  equipment  is  capable  of  rf  voltages  high  enough  so  that  u  >  uj 

5  p 

However,  when  such  voltages  are  applied,  the  observed  plasma  oscillations  prevent 
any  transit  time  measurements.  For  this  eason  most  experiments  have  been  per¬ 
formed  below  this  limit.  Fundamental  changes  are  observed  when  the  field  gradient 
approaches  this  limit.  These  changes  alter  the  nature  of  the  acceleration.  In  an 
effort  to  understand  these  changes  an  harmonic  analysis  that  includes  the  fundamental 
frequency  and  the  first  harmonic  but  neglects  higher  frequency  terms  has  been  per¬ 
formed.  The  results  provide  some  insight  into  changes  that  occur  as  the  limit  is 
approached,  but  are  not  valid  for  <u  >  <u  . 

o  r 

The  plasma  is  assumed  to  be  made  up  of  two  rigid  slabs  of  ions  and  electrons ,  and 
the  differential  equations  to  be  solved  are  * 

d28  d8  9  0  „  ,  , 

rrm  +  V  ““  +  6J  B  +  6J  _  B  COS  (co t  +  ~  —  I  E*  I  cos  (cot  +  y» )  (C“i) 

dt2  dt  P  g  m 
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where  the  plasma  is  assumed  to  be  in  '  localized  region  with  the  center  of  mass  at 
and  the  rf  field  magnitude  is  Iff/. 


According  to  Equation  C-2,  the  average  acceleration  of  the  plasma  during  one  period 


of  the  rf  is 


N  it2  '  H  *  \ 


(C-3) 


To  calculate  the  right  hand  side,  it  is  necessary  to  obtain  the  separation  5  .  if 

5  is  expressed  as  a  Fourier  series,  the  only  component  that  is  important  is  that 
of  fundamental  frequency,  since  the  other  components  will  disappear  when  the  aver¬ 
age  of  Equation  C-3  is  taken. 

The  separation  8  is  assumed  to  be  in  the  form 


5  =l|ff|  A  _J _ +  B  coscot  +  C  cos  (2  cot  +  <f>  +  {J,) 


( C-4 ) 


where  the  phase  oi  the  fundamental  frequency  is  included  in  the  phase  factor  in  the 
differential  equations.  Equation  C-4  is  now  substituted  into  Equation  C-i,  and  the 
collection  of  terms  leads  to  the  five  equations 
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sin  >ot; 
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while  the  higher  frequency  term  (fat )  is  dropped.  These  five  equations  can 
be  solved  for  the  five  quantities  A,  B,  C,  <f>,  and  0  , 


The  result  of  solving  Equation  (08)  and  (09)  for  0  and  B/C  is 

ton  0  =  —  - '•••Bwy — —  (C-10) 
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The  result  of  solving  equation  (06)  is: 
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The  above  quantities  are  substituted  into  Equations  (03)  and  (07),  the  result  being 
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Finally  substituting  back  into  Equation  (C-3)  the  average  plasma  acceleration  is: 


where 

*  (a>2  -  +  (toiv)*  (C-17) 

P 

If  (o4  is  very  small  compared  to  ^  ,  the  analysis  is  similar  to  the  origin,/ 
theory,  and  Equation  C-16  reduces  to  an  equation  similar  tr»  Equation  7.  Wu^». 
the  applied  frequency  is  close  to  the  critical  frequency  and  ine  electron  collision 
frequency  is  small,  the  limit  of  applicability  of  the  original  theory  is 

S.  «  S  -  J  (C-18) 

dZ  «  p 

For  an  average  field  gradient  of  dE/dz  =  KE  where  K  e  1M  for  the 

present  apparatus,  it  is  possible  to  calculate  the  theoretical  final  plasma  velocity 
from  Equation  (C-16).  One  curve  is  shown  in  Figure  19  and  is  compared  to  the 
corresponding  curve  from  the  simple  theory.  For  this  particular  case  a  deviation 
of  10%  in  the  plasma  velocity  occurs  when  the  field  gradient  is  only  0. 3  of  m/e  . 

The  analysis  shows  that  the  effects  of  large  field  gradients  are  important  for  the 
experiments.  It  does  not  help  in  explaining  the  differences  between  theory  and 
experiment  shown  in  Figure  16,  since  the  improved  theory  tends  to  increase  the 
discrepancy,  rather  than  reduce  it. 

It  should  be  noted  that  the  application  of  this  analysis  has  been  made  only  for  the 
region  beyond  the  field  maximum,  i.e.  ,  the  main  accelerating  region.  It  has  been 
assumed  that  the  transition  between  the  zero  field  inside  the  plasma  source  and  the 
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VELOCITY  (I03cm/SEC) 


Figure  19.  Effect  of  Large  Field  Gradients  on  Plasma  Velocity 
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maximum  field  outside  the  orifice  is  so  sharp  that  there  are  no  acceleration  effects. 
This  is  equivalent  to  assuming  that  the  field  gradients  in  tnis  region  are  very  large. 
The  experiments  discussed  in  Appendix  E  indicate  that  this  assumption  fails  to  hold 
for  rf  voltages  below  iOO  volts.  At  these  low  voltages  deceleration  may  take  place 
before  the  plasma  enters  the  maximum  field  region,  and  thus  might  decrease  the 
final  plasma  velocities  observed.  Such  aneffectwould  increase  the  observed  slope 
of  plasma  velocity  vs.  rf  field,  and  might  explain  the  slope  discrepancy  between 
the  present  theory  and  experiment  which  is  shown  in  Figure  16. 
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APPENDIX  D 


ION  TEMPERATURE  DETERMINATION  FROM  AMBIPOLAR  DIFFUSION 


Plasma  generated  by  the  mercury  Igniter  source  moves  into  the  rf  interaction 
region  by  ambipolar  diffusion.  A  measurement  of  the  ambipolar  diffusion  velocity 
together  with  the  plasma  electron  temperature  will  give  a  determination  of  the  plas¬ 
ma  ion  temperature. 


The  plasma  diffusion  velocity,  «  ,  fora  two-charge  system  is 


19 


v*K~  +  tf  K+ 

V  m  ■  1  — 

K~ 


(D-l) 


where  (+)  denotes  ions  and  (-)  denotes  electrons,  «  is  the  diffusion  velocity  of 
particles  in  the  absence  of  oppositely  charged  particles,  and  K  is  the  mobility. 
Since  the  mobility  of  the  ions  is  much  less  than  that  of  the  electrons,  the  diffusion 
velocity  is  approximately 


V  m 


<D-2) 


where  denotes  thermal  velocity. 

20 

By  the  Einstein  relationship, 

K+  D+T~ 

K~  D~T+ 


(D-3) 


where  D  is  the  diffusion  constant  of  the  particles  and  T  is  the  temperature.  The 
diffusion  constant  can  be  expressed  in  terms  of  the  mean  free  path  L  and  the 
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thermal  velocity 
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The  mean  free  path  for  Hg  lone  in  a  Hg  gaa  ia 


(D-4) 
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VT  Lg 

and  for  20, 000°K  electroni  in  a  Hg  gaa 
L~  •  2.8  Lg 

where  Lg  la  the  mean  free  path  of  the  neutral  Hg  atoma. 

The  ratio  of  the  ion  to  the  electron  diffuaion  oenatant  ia 

0+  VT  v+  v+ 

_  m  1  m  £  1 

O'  2.8  vT-  v- 

Comblning  equationa  (2),  (3)  and(7) 

fr) 

The  average  thermal  velocity  of  the  iona  la 

* J. ZT 

r  *  2  *  M 

where  ia  the  maaa  of  the  iona  4  la  the  Boltsman  oonatant. 


(D-5) 


(D-«) 


(D-7) 


(D-8) 


(D-fl) 


The  amblpolar  diffuaion  velooity  la  therefore,  a  funotion  of  the  electron  and  ion 
temperatures  and  the  mass  of  the  iona, 


&('-  fj 


(D-10) 
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Solving  for  T + 


When 
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APPENDIX  E 


EFFECT  OF  RF  VOLTAGE  ON  MAGNITUDE  OF  PROBE  CURRENT 

A.  INTRODUCTION 

Experiments  on  probe  currents  at  low  rf  voltages  yielded  the  following  re¬ 
sults.  The  peak  current  on  a  probe  decreased  (by  a  factor  of  ten)  as  the  rf  volt¬ 
age  was  increased,  and  after  reaching  a  minimum  then  increased  with  voltage. 

The  experimental  conditions  are  describt  in  section  IV  and  the  results  are  shown 
in  Figure  15.  These  experimental  results  can  be  qualitatively  explained  as  follows; 
As  the  rf  voltage  increases,  the  field  barrier  (  proportional  to  the  square  of  the 
maximum  field)  increases  in  the  region  outside  the  source  orifice.  The  number 
of  electrons  that  have  sufficient  thermal  energy  to  cross  the  barrier  decreases  as 
the  barrier  is  raised.  For  higher  voltages,  electrons  of  the  proper  phase  can  be 
accelerated  all  the  way  to  the  maximum  field  point  in  less  than  one  rf  cycle,  and 
are  not  contained  by  the  field  barrier;  hence,  the  number  of  electrons  observed 
increases  for  higher  rf  voltages.  Not  only  the  decrease  and  increase  of  electrons 
with  rf  voltage,  but  also  the  magnitude  of  the  rf  voltage  at  which  the  minimum  is 
observed,  can  be  explained  by  the  theory. 

The  decrease  of  intensity  with  increasing  rf  voltage  and  the  later  increase  are  due 
to  different  effect,  and  the  observed  sharp  minimum  is  partly  a  coincidence.  Under 
different  circumstances  the  minimum  might  be  quite  broad,  and  the  electrons  might 
be  contained  by  the  field  barrier  for  several  orders  of  magnitude  change  in  rf  volt¬ 
age.  The  theory  is  therefore  divided  in  two  parts,  the  first  explaining  the  decrease, 
and  the  second  the  increase. 
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The  theory  Is  based  on  the  behavior  of  the  electrons,  and  neglects  their  interactions 
with  the  ions .  A  theory  that  includes  the  electron-ion  interactions  in  the  region  in 
which  the  plasma  enters  the  rf  field  has  not  yet  been  formulated,  but  there  is  reason 
to  believe  that  for  the  data  shown  in  Figure  15  these  interactions  are  not  important. 


B.  DECREASE  IN  INTENSITY 

The  deceleration  of  electrons  for  an  applied  frequency  much  greater  than  the 
plasma  frequency,  is  given  by  theory 1  as 

fi  . .  jL  JL  i£  <E-i) 

di*  2m2  ai2  dZ 

If  an  electron  starts  inside  the  Igniter  source  with  an  initial  velocity  vi  in  a  region 
of  no  electric  field,  Integration  of  equation  E-l  shows  that  this  velocity  must  be 
greater  than  a  certain  value,  depending  on  the  maximum  electric  field  Emax  ,  in 
order  that  the  electrons  pass  the  field  maximum. 


>/7  m 


max 
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(E-2) 


It  is  assumed  that  the  electron  velocities  have  a  Maxwell  distribution, 
for  velocities  in  one  specified  direction  is 


2KT 


m 


,  which 


(E-3) 


where  T  is  the  temperature  of  the  electrons .  The  observed  density  of  electrons 
will  be  proportional  to  the  number  of  electrons  with  a  velocity  greater  than  the  min¬ 


imum,  i.e. , 


00  M 


(E-4) 
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where 


(E-5) 


If  tho  density  of  electrons  N  is  plotted  as  a  function  of  x  ,  the  result  as  shown  in 
Figure  20  is  the  familiar  integral  of  the  error  curve.  If  the  slope  at  a?  -  0  is 
extrapolated  to  the  x-axis,  the  Intercept  is  a  characteristic  %e  ,  which  is  equal  to 


“  0)o  n/2  \frT 

(dl/dx)^  n//fT  2 


(E-6) 


The  curve  for  N  vs.  maximum  electric  field  would  have  exactly  the  same  shape, 
with  the  scale  factor  for  the  x-axis  being  given  in  Equation  E-5.  Solving  for  the 
characteristic  electric  field,  the  result  is 


where  and  e  are  the  mass  and  charge  of  the  electron  and  k  is  the  Boltzmann 
constant.  In  the  present  experiment  the  applied  frequency  was  w/2  n  =  HO  Me  ,  and 
the  measured  electron  temperature  T  was  20,000°K  as  reported  in  section  IV.  The 
calculated  characteristic  field  is  then: 


Ee  -  47  v/cm  (E-8) 

This  field  is  equivalent  to  ar.  rf  potential  of  55  volts  in  the  present  experiment.  If 
the  first  part  of  the  curve  in  Figure  15  is  extrapolated,  the  x-axis  intercept  is  about 
60  volts,  which  is  in  agreement  with  the  theory. 

C.  INCREASE  IN  INTENSITY 

The  theory  of  the  electric  field  barrier  assumes  that  the  particle  makes  several 
oscillations  which  result  in  a  net  deceleration.  If  the  amplitude  of  oscillation  be¬ 
comes  of  the  same  order  of  magnitude  as  the  extent  of  the  region  of  the  gradient, 
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electrons  that  arrive  during  a  certain  phase  of  the  rf  cycle  will  not  be  decelerated. 
The  amplitude  of  the  electron  oscillations,  if  ion  interactions  are  neglected,  is 

1 5  I  -  —  —  u»oj  (E-9) 

m  CO2 

The  distance  D,  in  which  the  electric  field  increases  from  a  small  value  such  as  1% 
to  its  maximum  value  outside  the  orifice,  is  0. 1  cm,  as  can  be  seen  in  Figure  5. 
Setting  this  distance  equal  to  |S|  above,  the  equation  can  be  solved  for  the  electric 
field  at  which  electrons  will  start  to  leak  through  the  potential  barrier. 

E  m  =  4Sv/cm  (E-10) 

e 

This  field  at  which  electrons  begin  to  leak  through  the  rf  barrier  is  equivalent  to  an 
rf  potential  of  55  volts.  This  is  an  excellent  agreement  with  Figure  15;  as  the  rf 
potential  is  increased  beyond  this  limit,  a  larger  number  of  electrons  pass  through 
the  barrier. 


D.  CONCLUSIONS 


Equation  E-10  should  be  compared  with  Equation  E-7. 


Decrease: 


47  v/ cm 


Increase: 
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Dmo)2 
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45  v/cm 


(E-7) 


(E-10) 


If  any  one  of  the  factors  in  the  equation  i3  changed,  the  decrease  and  increase  would 
no  longer  be  at  the  same  field,  and  the  shape  of  the  intensity  vs.  voltage  curve  would 
be  quite  different.  If  Equation  E-10  had  had  a  much  smaller  numerical  result  than 
Equation  E-7  the  probe  current  would  not  drop  by  a  factor  of  five,  and  if  it  had  been 
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higher,  the  minimum  would  be  much  broader.  The  minimum  should  inde  d  be¬ 
come  broader  if  the  orifice  diameter  or  the  applied  frequency  is  increased,  while 
it  will  be  reduced  in  magnitude  if  the  initial  electron  temperature  is  increased. 

The  theory  can  be  compared  with  the  experimental  data  shown  in  Figure  15.  The 
theory  is  in  good  agreement  as  to  the  sharp  minimum  observed,  and  the  value  of 
the  electric  field  at  which  this  takes  place.  There  is  also  reasonable  agreement 
between  the  shape  of  the  curve  in  Figure  20  and  the  decrease  in  intensity  shown 
in  Figure  15. 
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